skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Stockli, Daniel F"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Halogens (F, Cl, Br, I) are primary components of volcanic gas emissions and play an essential role in continental arc magmatic environments due to their solubility in fluids that generate metallic ore deposits. Despite their ubiquity, the behavior and budget of halogens in continental arc environments are poorly constrained. We investigated the plutonic and volcanic halogen budgets in intermediate-to-felsic igneous rocks (56–77 wt% SiO2) from the Sierra Nevada (California) - a Mesozoic continental arc where plutonic and volcanic outcrops can be correlated via their geographic, compositional, and geochronologic framework. We measured the halogen concentrations of bulk rock powders and their leachates via ion chromatography (F, Cl) and ICP-MS (Br, I). Halogen concentrations in our rock powders range between 107–727 μg/g F, 13–316 μg/g Cl, 2–323 ng/g Br, and 1–69 ng/g I. In contrast, leachates yielded 3–4 orders of magnitude less Cl and F, one order of magnitude less I, and similar amounts of Br compared to their corresponding bulk rocks. Preliminary data show no significant differences between volcanic and plutonic samples, suggesting that halogen concentrations in these rocks are insensitive to shallow fractionation. Although F and I exhibit no correlation with major element compositions, Cl and Br display negative trends with increasing SiO2 and K2O, and positive trends with increasing Fe2O3T, MnO, MgO, CaO, and TiO2, suggesting mafic minerals as important hosts of structurally bound halogens. Overall, Sierran plutonic rocks display low halogen contents (max. F, Cl = 727, 315 μg/g), consistent with biotite- and apatite-bearing granitoids reported in [1]. This work suggests that halogens do not preferentially enrich in shallow plutonic or volcanic portions of a continental arc system and that mafic mineral phases likely serve as primary reservoirs of these elements in intermediate-to-felsic igneous rocks. These hypotheses will be further investigated in future work through in-situ analysis of halogen concentrations in crystals. [1] Teiber, Marks, Wenzel, Siebel, Altherr & Markl (2014), Chemical Geology, vol. 374–375, pp. 92–109, doi: 10.1016/j.chemgeo.2014.03.006 
    more » « less
  2. Sedimentary basins record crustal-scale tectonic processes related to the construction and demise of orogenic belts, making them an invaluable archive for the reconstruction of the evolution of the North American Cordillera. In southwest Montana, USA, the Renova Formation, considered to locally represent the earliest accumulation following Mesozoic−Cenozoic compressional deformation, is widespread but remains poorly dated, and its origin is debated. Herein, we employed detrital zircon U-Pb and (U-Th)/He double dating and sanidine 40Ar/39Ar geochronology in the context of decimeter-scale measured stratigraphic sections in the Renova Formation of the Muddy Creek Basin to determine basin evolution and sediment provenance and place the basin-scale record within a regional context to illuminate the lithospheric processes driving extension and subsidence. The Muddy Creek Basin is an extensional half graben in southwest Montana that is ∼22 km long and ∼7 km wide, with a >800-m-thick sedimentary package. Basin deposition began ca. 49 Ma, as marked by multiple ignimbrites sourced from the Challis volcanic field, which are overlain by a tuffaceous fluvial section. Fluvial strata are capped by a 46.8 Ma Challis ignimbrite constrained by sanidine 40Ar/39Ar dating. An overlying fossiliferous limestone records the first instance of basinal ponding, which was coeval with the cessation of delivery of Challis volcanics−derived sediment into the Green River Basin. We attribute initial ponding to regional drainage reorganization and damning of the paleo−Idaho River due to uplift and doming of the southern Absaroka volcanic province, resulting in its diversion away from the Green River Basin and backfilling of the Lemhi Pass paleovalley. Detrital zircon maximum depositional ages and sanidine 40Ar/39Ar ages show alternating fluvial sandstone and lacustrine mudstone deposition from 46 Ma to 40 Ma in the Muddy Creek Basin. Sediment provenance was dominated by regionally sourced, Challis volcanics−aged and Idaho Batholith−aged grains, while detrital zircon (U-Th)/He (ZHe) data are dominated by Eocene cooling ages. Basin deposition became fully lacustrine by ca. 40 Ma, based on an increasing frequency of organic-rich mudstone with rare interbedded sandstone. Coarse-grained lithofacies became prominent again starting ca. 37 Ma, coeval with a major shift in sediment provenance due to extension and local footwall unroofing. Detrital zircon U-Pb and corresponding ZHe ages from the upper part of the section are predominantly Paleozoic in age, sourced from the Paleozoic sedimentary strata exposed in the eastern footwall of the Muddy Creek detachment fault. Paleocurrents shift from south- to west-directed trends, supporting the shift to local sources, consistent with initiation of the Muddy Creek detachment fault. Detrital zircon maximum depositional ages from the youngest strata in the basin suggest deposition continuing until at least 36 Ma. These data show that extension in the Muddy Creek Basin, which we attribute to continued lithospheric thermal weakening, initiated ∼10 m.y. later than in the Anaconda and Bitterroot metamorphic core complexes. This points to potentially different drivers of extension in western Montana and fits previously proposed models of a regional southward sweep of extension related to Farallon slab removal. 
    more » « less
  3. The uplift history of the Sierra Nevada, California, is a topic of long-standing disagreement with much of it centered on the timing and nature of slip along the range-bounding normal fault along the east flank of the southern Sierra Nevada. The history of normal fault slip is important for characterizing the uplift history of the Sierra Nevada, as well as for characterizing the geologic and geodynamic factors that drove, and continue to drive, normal faulting. To address these issues, we completed new structural studies and extensive apatite (U-Th)/He (AHe) thermochronometry on samples collected from three vertical transects in the footwall to the east-dipping southern Sierra Nevada normal fault (SNNF). Our structural studies on bedrock fault planes show that the SNNF is a steeply (~70°) east-dipping normal fault. The new AHe data reveal two elevation-invariant AHe age arrays, indicative of two distinct periods of cooling and exhumation, which we interpret as initiation of normal faulting along the SNNF at ca. 28–27 Ma with a second phase of normal faulting at ca. 17–13 Ma. We argue that beginning in the late Oligocene, the SNNF marked the now long-standing stable western limit, or break-away zone, of the Basin and Range. Slip along SNNF, and the associated unloading of the footwall, likely resulted in two periods of uplift of Sierra Nevada during the late Cenozoic. Trench retreat, driven by westward motion of the North American plate, along the Farallon–North American subduction zone boundary, as well as the gravitationally unstable northern and southern Basin and Range pushing on the cold Sierra Nevada, likely drove the late Oligocene- aged normal slip along the SNNF and the similar-aged but generally local and minor extension within the Basin and Range. We posit that the thick proto–Basin and Range lithosphere was primed for late Oligocene extension by replacement of the steepening Farallon slab with hot and buoyant asthenosphere. While steepening of the Farallon slab had not yet reached the southern Sierra Nevada by late Oligocene time, we speculate that late Oligocene slip along the SNNF reactivated a late Cretaceous dextral shear zone as the Sierra Nevada block was pulled and pushed westward in response to trench retreat and gravitational potential energy. The dominant middle Miocene normal fault-slip history along the SNNF is contemporaneous with high-magnitude slip recorded along range-bounding normal faults across the Basin and Range, including the east-adjacent Inyo and White mountains, indicating that this period of extension was a major regional tectonic event. We infer that a combination of slab-driven trench retreat along the Juan de Fuca–North America subduction zone boundary and clockwise rotation of the southern ancestral Cascade Range superimposed on continental lithosphere pre-conditioned for extension drove this episode of middle Miocene normal slip along the SNNF and extension to the east across the Basin and Range. Transtensional plate motion along the Pacific–North America plate boundary, and likely a growing slab window, continued to drive extension along the SNNF and the western Basin and Range, but not until ca. 11 Ma when the Mendocino triple junction reached the latitude of our northernmost (U-Th)/He transect. 
    more » « less
  4. Evidence from landscape evolution may provide critical constraints for past geodynamic processes, but has been limited by the large uncertainties of topographic reconstructions. Here we present continuous 30-million-year rock uplift histories for three catchments in the Calabrian forearc of southern Italy, using a data-driven inversion of tectonic geomorphology measurements. We find that rock uplift rates were high (>1 mm yr−1) from about 30 to 25 million years ago (Ma) and progressively declined to <0.4 mm yr−1 by ~15 Ma, then remained low before abruptly increasing around 1.5–1.0 Ma. These uplift rates do not match the forearc’s subduction velocity record, implying that uplift was not dominated by crustal thickening due to subduction-driven sediment influx. Through comparisons with slab descent reconstructions, we instead argue that the forearc uplift history primarily reflects the progressive establishment and abrupt destruction of an upper-mantle convection cell with strong negative buoyancy. We suggest that the convection cell vigour increased as the slab-induced mantle flow field began to interact with the 660-km mantle transition zone, causing uplift rates to decline from 25 to 15 Ma. Then, once the slab encountered the transition zone, the fully established convection cell subdued uplift rates, before being disrupted by slab fragmentation in the Quaternary, driving rapid forearc uplift. 
    more » « less
  5. Crustal thickening along the Indus-Yarlung suture zone during the India-Asia collision was primarily accommodated by slip along the north-dipping Gangdese thrust (ca. 27–18 Ma) and south-dipping Great Counter thrust (ca. 25–10 Ma). However, the along-strike continuities, geometries, and timings of these thrusts remain unclear, resulting in an inadequate understanding of Himalayan-Tibetan orogenesis. In this study, we performed geologic mapping, strain analyses, geo/thermochronology, and thermobarometry across the easternmost Himalayan orogen (i.e., the northern Indo-Burma thrust belt), specifically: (1) the Tidding thrust and easternmost Indus-Yarlung suture zone (i.e., Tidding mélange complex) in its hanging wall; and (2) the Lohit thrust and Jurassic–Cretaceous Gangdese batholith and Mesoproterozoic basement (i.e., Lohit Plutonic Complex) in its hanging wall. The Tidding thrust is a north-dipping, top-south mylonitic shear zone that was active by ca. 36–30 Ma, during which hanging-wall mélange rocks were exhumed from ~33–38 km depth. The geometry, kinematics, and initiation age of the Tidding thrust contrast those of the top-north Great Counter thrust at the same structural position to the west. North of the Tidding thrust, the Lohit thrust is a ~5-km-wide, subvertical, north-side-up mylonitic shear zone that contains a basal, discrete “Lohit thrust fault”. Results of electron backscatter diffraction analyses across the Lohit thrust shear zone show that deformation fabric intensity and finite strain magnitudes decrease southwards toward the discrete thrust fault. This spatial relationship may be the result of transient peak strain during the lifespan of the shear zone. The Lohit thrust was active by ca. 25–23 Ma, during which hanging-wall basement and batholithic root rocks were exhumed to mid-crustal depths. The Lohit thrust and Gangdese thrust to the west are located at the same structural position and have comparable geometries, kinematics, and timings. Based on these similarities and previous findings, we interpret that the Lohit and Gangdese thrusts are correlative segments of a single, orogen-wide thrust system that accommodated crustal thickening along the Indus-Yarlung suture zone during the Oligocene–Miocene. 
    more » « less
  6. Crustal thickening has been a key process of collision-induced Cenozoic deformation along the Indus-Yarlung suture zone, yet the timing, geometric relationships, and along-strike continuities of major thrusts, such as the Great Counter thrust and Gangdese thrust, remain inadequately understood. In this study, we present findings of geologic mapping and thermo- and geochronologic, geochemical, microstructural, and geothermobarometric analyses from the easternmost Indus-Yarlung suture zone exposed in the northern Indo-Burma Ranges. Specifically, we investigate the Lohit and Tidding thrust shear zones and their respective hanging wall rocks of the Lohit Plutonic Complex and Tidding and Mayodia mélange complexes. Field observations are consistent with ductile deformation concentrated along the top-to-the-south Tidding thrust shear zone, which is in contrast to the top-to-the-north Great Counter thrust at the same structural position to the west. Upper amphibolite-facies metamorphism of mélange rocks at ∼9−10 kbar (∼34−39 km) occurred prior to ca. 36−30 Ma exhumation during slip along the Tidding thrust shear zone. To the north, the ∼5-km-wide Lohit thrust shear zone has a subvertical geometry and north-side-up kinematics. Cretaceous arc granitoids of the Lohit Plutonic Complex were emplaced at ∼32−40 km depth in crust estimated to be ∼38−52 km thick at that time. These rocks cooled from ca. 25 Ma to 10 Ma due to slip along the Lohit thrust shear zone. We demonstrate that the Lohit thrust shear zone, Gangdese thrust, and Yarlung-Tsangpo Canyon thrust have comparable hanging wall and footwall rocks, structural geometries, kinematics, and timing. Based on these similarities, we interpret that these thrusts formed segments of a laterally continuous thrust system, which served as the preeminent crustal thickening structure along the Neotethys-southern Lhasa terrane margin and exhumed Gangdese lower arc crust in Oligocene−Miocene time. 
    more » « less